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Abstract 



In this paper, we study nonleptonic charmless B decays to two hght 
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the contributions from the chirally enhanced power corrections and 
weak annihilation. Predictions for the CP-averaged branching ratios 
and CP-violating asymmetries are given. Within the reasonable range 
of the parameters, we find that our predictions for the branching ra- 
tios of -B ^ PP are consistent with the present experimental data. 
But because of the logarithmic divergences at the endpoints in the 
hard spectator scatterings and weak annihilation, there are still large 
uncertainties in these predictions. 
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I. INTRODUCTION 



In the standard model (SM), CP violation and quark mixing are closely related to 
each other. The unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) matrix implies 
various relations among its elements. The most commonly studied is the "unitarity tri- 
angle": VudV*i^ + Vcd^*^ + VtdVtl = 0. The study of B meson decays is mainly to make 
enough independent measurements for the sides and angles (a, P and 7) of this unitarity 
triangle. For example, in principle, sin 2/? can be determined by measurements of the 
time dependent CP asymmetries of the decay —>■ J/'^Ks and other related modes 
with an uncertainty less than 1%. On the other hand, we can also extract or give a 
constraint for the angle a from the decay B — »• tttt, and the angle 7 from B tctt and 
TcK. Compared with determination of sin 2(3, the extractions of a and 7 from exclusively 
non-leptonic charmless B decays suffer from significant theoretical uncertainties for the 
strong interactions. 

Experimentally, many B experiment projects have been running (CLEO, BaBar, Belle, 
etc.), or will run in forthcoming years (BTeV, CERN LHCb, DESY HeraB, etc.). Several 
collaborations have reported their latest results recently H]-^, and more B decay channels 
will be measured with great precision soon. Among these decay channels, quite a few 
are nonleptonic charmless two-body modes, such as -B ^ nK, tttt, Krj', etc. With the 
accumulation of the experimental data, the theorists are urged to gain a deep insight into 
rare hadronic B meson decays, and to reduce the theoretical errors in determining the 
CKM parameters from the experimental data. 

Theoretically, the nonleptonic charmless two-body B decays have been widely dis- 
cussed and carefully studied in the framework of the effective Hamiltonian and factor- 
ization hypothesis. The effective Hamiltonian is obtained by operator product expansion 
(OPE) and the renormalization group (RG) method, and generally expressed as the prod- 
uct of the Wilson coefficients and the effective operators, in forms of the four-quarks and 
magnetic moment operators. The Wilson coefficients can be calculated reliably by the 
perturbation theory. In the SM, they have been evaluated to next-to- leading order 
Thus, the main task for us is to compute the hadronic matrix elements of the effective 
operators. However, the complexity of QCD dynamics does not allow us to compute them 
directly from first principles. Generally, we resort to the factorization approximation ^ , 
ignoring the soft interactions between the ejected meson and recoiling system. Under this 
approximation, the hadronic matrix element of the effective operator can be parameter- 
ized into the product of the meson decay constant and the meson-meson transition form 
factor. A comprehensive analysis for nonleptonic charmless two-body B decays based on 
the naive factorization and/or generalized factorization was done in the past decade, and 
achieved great success. 

The factorization approximation does hold in the limit that the soft interactions in the 
initial and final states can be ignored. It seems that the argument of color-transparency 
can give reasonable support to the above limit . Because the b quark is heavy, the quarks 
from b quark decay move so fast that a pair of quark and antiquark in a small color-singlet 
object decouples from the soft interactions. But for the case of naive factorization the 
shortcomings are obvious. First, the renormalization scheme and scale dependence in the 
hadronic matrix elements of the effective operators are apparently missed. Then the full 
decay amplitude predicted by the Bauer-Stech-Wirbel (BSW) model remains dependent 
on the renormalization scheme and scale, which are mainly from the Wilson coefficients. In 
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recent years, many researchers have generahzed the naive factorization scheme and made 
remarkable progress, such as the scheme and scale independent effective Wilson coeffi- 
cients, effective color number, which is introduced to compensate the "non-factorizable" 
contributions, etc. |^. Furthermore, some progress in nonperturbative methods, such as 



lattice QCD, QCD sum rules etc. ||8|-[10|, allow us to compute many non-perturbative pa- 
rameters in B decays, such as the meson decay constants and meson-meson transition form 
factors. Every improvement allows us to have a closer look at the B nonleptonic decays. 
In []Tll , p!2| , the authors gave a comprehensive and detailed analysis for the nonleptonic 
two-body charmless B decays within the frame of the generalized factorization. 

Two year ago, Beneke, Buchalla, Neubert, and Sachrajda (BENS) gave a QCD factor- 
ization formula in the heavy quark limit for the decays 5 — >■ vrvr |]13[. They pointed out 
that the radiative corrections from hard gluon exchange can be systematically calculated 
by use of the perturbative QCD method in the heavy quark limit, i.e., neglecting the 
power contributions of KqcD/^b- This factorization formula can be justified in the case 
that the ejected meson from the h quark decay is a light meson or an onium, no matter 
whether the other recoiling meson which absorbs the spectator quark in the B meson is 
light or heavy. But for the case that the ejected meson is in an extremely asymmetric con- 
figuration, such as a D meson, this factorization formula does not hold. The contributions 
from hard scattering with the spectator quark in the B meson are also involved in their 
formula. This kind of contribution cannot be contained in the naive factorization, but it 
appears at the order of a^. So the naive factorization can be recovered if one neglects the 
radiative corrections and suppressed power KQco/f^b contributions in the QCD factor- 
ization. The non-factorizable contributions in the naive factorization can be calculated 
perturbatively, so we do not need a phenomenological parameter N^^^ to compensate the 
non-factorizable effects. Moreover, because the h quark mass is not very large in reality, 
the potential power corrections, in particular, the chirally enhanced corrections, will play 
an important role in i? — > PP |I4|-[T6|]. So the twist-3 light-cone wave function of the 
light pseudoscalar must be taken into account. Unfortunately, QCD factorization breaks 
down at twist-3 level. There is a logarithmic divergence in the hard spectator scattering 
at the endpoint of the twist-3 light-cone distribution amplitude (LCDA). 

Phenomenologically, the QCD factorization (BBNS approach) has been applied to 
study many B meson decay modes, such as B ^ D^*^- [|17|,|18l, vrvr, tiK lT||gjl|Jl|], 
and other interesting channels pTl-121]. In particular, in a recent paper |TB[, Beneke et 



al. gave a detailed phenomenological analysis for B ^ irir and nK within the frame 
of QCD factorization, in which the contributions from weak annihilation are also taken 
into account. In principle, weak annihilation is power suppressed within the frame of 
QCD factorization. Nevertheless, as emphasized in Ref. these contributions may be 
numerically important for the decays B — > PP. An explicit calculation shows that there 
are very significant endpoint divergences in weak annihilation. This indicates that the 



weak annihilation may be dominated by the soft contributions. In Ref. |T^, the authors 
give a phenomenological treatment of the endpoint divergences in both the hard spectator 
scattering and the weak annihilations. Such a treatment is really a model- dependent 
procedure. However, it may be helpful to estimate the importance of annihilation. The 
numerical results show that weak annihilation indeed gives evident corrections to the 



branching ratios of the decays B ^ titt and ttK in Ref. |jT^, and also causes very large 



uncertainties in the prediction. Assuming the universality of the endpoint divergence, the 
authors of Ref. [|16| also gave a global fit for the CKM parameters by use of the branching 
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ratios of i? — tttt and tvK in the frame of QCD factorization, and derived a new constraint 
for the (p, fj) complex plane. The new constraint is slightly different from that derived 
by the other standard fit, but they are compatible with each other. 

In this work, we will apply QCD factorization to all modes of B decay into two light 
pseudoscalar mesons, and give a careful and detailed analysis. For the decay modes 
B ^ nir and ttK, our work is actually a re-check of Beneke et al. , and we find that 
our results are consistent with theirs. The B — >■ P?]^'-* has also been studied by M. Z. Yang 



and Y. D. Yang within the frame of QCD factorization The differences between their 



work and ours are that the other type twist-3 light-cone distribution amplitude (j)„ is taken 
into account, which guarantees gauge invariance of our results; and the contributions from 
weak annihilation are also considered in our work. In addition, the digluon mechanism 
for the production of r^^'^ is taken into account in Ref. I^Tl- Generally, it is believed that 
the anomalous coupling of T]'g*g* can account for the large branching ratio of B ^ Krj' . 
However, whether the form factor i]^'^g*g* is pertubatively calculable or not is still an 
open question. In particular, at the endpoint of the light-cone distribution amplitudes 
(LCDAs), where the virtualities of the two gluons are both soft, the validity of the sup- 
pression from the LCDAs is questionable. Therefore, in our work, digluon mechanism is 
not considered although it can enhance the branching ratio of i? — > Krj'. For B KK, 
our work is totally new; the study of these decay modes is expected to give a constraint 
to the endpoint divergence in weak annihilation. 

This paper is organized as follows. In Sect. II, we give a brief review of the effective 
Hamiltonian and QCD factorization including the chirally enhanced power corrections 
and weak annihilations. For a quantitative analysis of -B ^ PP, the necessary input 
parameters are discussed in Sect. III. In Sect. IV, the numerical results for the CP- 
averaged branching ratios and CP asymmetries for B PP are given and Sect. V is 
devoted to the conclusions. 



II. THEORETICAL FRAME FOR B RARE DECAYS 
A. The eflfective Hamiltonian 

B decays involve three characteristic scales which are strongly ordered: m\Y ^ S> 
Aqcd- How to separate or factorize these three scales is the most essential question in 
studying B hadronic decays. 

With the operator product expansion method, the relevant \AB\ = 1 effective Hamil- 
tonian is given by M 



G p 

71 



/ 10 N 

q=u,c \ k=3 i 



+ h.c., (1) 



where Vq = V^bl/'^(for b d transition) or Vq = V^fel/*(for b s transition) and Cj(p) 
are the Wilson coefficients which have been evaluated to next-to-leading order with the 
perturbation theory and renormalization group method. 
In the Eq.(l), the four-quark operators Qi are given by 
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Ql = {Uaba)v-A{qi3Uf3)v-A 
Q2 = {uJ)p)v-A{.qpUa)v-A 
= {.qaK)v-All{.qf3(lp)v-A 



Q 

Q5 = (qak 



)V-A 



Eiq^q'/s] 



Q\ = iCaba)v-A{q/3C/3)v-A 
Q2 = {Cab/3)v-A{qi3Ca)v-A 

Qa = {qpba)v-AY.{qJi3)v-A 

g' 



V+A 



Q6 = 



V~Aj2iqaqf3)v+A 
q' 

Eeg'{qjp)v+A 



'a)v~A . 



Qy = l{qaba)v-AJ2eq'{qi3qi3)v+A Qs = |i 

q' q' 

Q9 = l{qaba)v-AT, eg'{qpq'p)v-A Qio = l{qpba)v-AT,eq'{q^q,3)v-A 

q' q' 



(2) 



and 

g g 

= -^^bqa(^'"' {I + l5)baF^u, QsG = ^mbq^a^'' t^^pil + -i-,)bpGl^, {q = d or s). 

(3) 

Here Ql and Ql are the tree operators, Q3 — Qq the QCD penguin operators, Qj — Qio 
the electroweak penguin operators, and Qy^ and Qsc the magnetic-penguin operators. 

In this effective Hamihonian for B decays, the contributions from the large virtual 
momenta of the loop corrections from the scale fi = 0{mb) to mw are given by the 
Wilson coefficients, and the low energy contributions are fully incorporated into the matrix 
elements of the operators Q. So the derivation of the effective Hamiltonian can be called 
"the first step factorization". 

Several years ago, the perturbative corrections to the Wilson coefficients in the SM 
were evaluated to next-to-leading order with the renormalization group method. We list 
the numerical results in the naive dimensional regularization (NDR) scheme and at the 
three scales mh/2, rrif, and 2mh in Table |. 



B. QCD factorization in the heavy quark limit 

After "first step factorization", the decay amplitude for B P1P2 can be written as 
A{B ^ P1P2) = Y.^iC,{^l){P^P2\Q^{^l)\B) . (4) 

i 

The central task now is to reliably evaluate the transition matrix element {PiP2\Qi{ij)\B) . 
Note that the matrix element contains two scales: 0{mb) and Aqcd- It is hoped that 
{PiP2\Qi{fi)\B) can be separated into the short- distance contributions relating to large 
scale 0{mh) and the long-distance contributions relating to the fundamental scale of 
QCD-Aqcd- Then the short-distance contributions, which are perturbatively calculable, 
should recover the scheme and scale dependence of the hadronic matrix elements, and the 
long-distance contributions can be parameterized into some universal non-perturbative 
parameters. 

In Ref. [0], Beneke, Buchalla, Neubert and Sachrajda show that this idea of factor- 
ization can be realized, at least at one- loop order, in the heavy quark limit. They argue 
that the emitted meson Pi carries large energy and momentum (about and there- 

fore can be described by leading-twist light-cone distribution amplitudes. It is natural 
to imagine that the soft gluons would decouple from the emitted Pi at leading order of 
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KqcD/fT^b since the qq pair in Pi form a small-size color dipole. However, the factor- 
ization requires that not only soft divergences, but also coUinear divergences, should be 
cancelled. Fortunately, explicit calculations at one-loop order show that, in the heavy 
quark limit, all infrared divergences vanish after summing over the four vertex correction 
diagrams [Figs. l(a)-l(d)]. As for the recoiled meson P2, because the spectator quark in 
the B meson is transferred to it as a soft parton, Beneke et al. believe that, instead of 
light-cone distribution amplitudes, one should use nonperturbative form factors F^^'^ to 
describe it. Q It should be noted that when the spectator quark interacts with one of 
the quarks in the emitted meson by a hard gluon exchange, the recoiling meson can be 
described by leading-twist distribution amplitudes. These hard spectator diagrams [Figs. 
1(g), 1(h)] can also be accounted for at the leading power of kQcolf^h- In summary, their 
factorization formula can be explicitly expressed as 

1 1 

{PiP^mB) = F''^''\^) J dxTl{x)^p,{x) + J didxdyTl\i,x,y)^B{i)^pAx)'^pM 



= (P1P2I Ji ® J2\B) ■ [1 + + 0{KQCD/m,)]. (5) 

We call this factorization formalism as QCD factorization or the BBNS approach. In 
Eq.(^, $b(0 and $Pi(a;) {i = 1,2) are the leading-twist wave functions of B and the 
light pseudoscalar mesons respectively, and T/'^^ denote hard-scattering kernels which 
are calculable in perturbation theory. At the order of a^, the hard kernels T^'^^ can be 
depicted by Fig. 1. Figures l(a)-l(d) represent vertex corrections. Figs. 1(e) and 1(f) 
penguin corrections, and Figs. 1(g) and 1(h) hard spectator scattering. 

Comparing this approach with naive factorization and/or generalized factorization, 
there are some interesting characteristics in the BBNS approach. 

(i) At the leading order of a^, it can reproduce the results of naive factorization; at 
the higher order of a^, the renormalization scheme and scale dependence for the hadronic 
matrix elements can be recovered from the hard-scattering kernels T/. The generalized 
factorization can also obtain the necessary dependence on scheme and scale. However, this 
dependence is based on one-loop calculations of quark-level matrix elements. According 
to Buras et al., quark-level matrix elements are accompanied by infrared divergences. 
To avoid such divergences, one usually assumes that external quark states are off shell. 
Unfortunately, this will introduce gauge dependence, which is also unphysical. In Ref. 



I4 , the authors give a calculation assuming the external quarks are on shell, and they 
show that the infrared singularities can be cancelled by assuming the final quarks are in 
the form of hadrons; their calculation is therefore gauge invariant. Their idea is almost the 



^Light-cone sum rules also justify this nonperturbative property for the B to light form factors. 
It is different from the conclusion of Keum et al. [25|. In Ref. |25|, the authors claim that the 
form factor F^'^{0) is perturbatively calculable in terms of the distribution amplitudes when the 
Sudakov JtT resummation and threshold resummation are taken into account. However, a recent 
detailed examination of such point of view by Descotes-Genon and Sachrajda shows that it is 
not justified [26|. They conclude that, the Sudakov form factors are not sufficient to suppress 
the contribution that comes from the nonperturbative region of large impact parameters, and 
therefore the form factor F^'^{0) is uncalculable in the standard hard-scattering approach. 
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same as the BBNS approach. The difference is that, in the BBNS approach, the ejected 
meson is represented as its hght-cone distribution amphtude; in Ref. [Q, the ejected 
qq pair is in such a configuration that it is a color-singlet object and the quark shares 
almost the same momentum as the anti-quark. Because results in Ref. [0 are still from 
a calculation at the quark level, some information from the distributions of momentum 
fraction in hadrons is ignored. 

(ii) Generalized factorization considers nonfactorizable contributions as intractable. 
Therefore, one may introduce one or more effective color numbers N^^^ to phenomeno- 
logically represent nonfactorizable contributions. Furthermore, N^^f is assumed to be 
universal to maintain predictive power. However, QCD factorization tells us that the 
nonfactorizable contribution is indeed factorizable and therefore calculable in the heavy 
quark limit. In consequence, N^^^ is calculable, and process dependent beyond the leading 
order of a^. 

(iii) An interesting result of QCD factorization is that strong phases come solely from 
hard scattering processes and are therefore calculable in the heavy quark limit. From 
Eq.(^, it is easy to conclude that the imaginary part of decay amplitude arises only from 
hard scattering kernels T/ because nonperturbative form factors and light-cone distribu- 
tion amplitudes are all real. T/ contains vertex corrections [Figs. 1 (a)-l(d)] and penguin 
corrections [Figs. 1(e), 1(f)]. Strong phases from penguin corrections are commonly called 
the Bander-Silverman-Soni (BSS) mechanism which is the unique source of strong 
phases in generalized factorization. However, the gluon virtuality of Fig. 1(e), which is 
well defined in QCD factorization, is ambiguous in generalized factorization and usually 
treated as a free parameter. In addition, vertex corrections will also contribute to strong 
phases in QCD factorization, which is missed in the generalized factorization. 

(iv) Hard spectator contributions [Figs. 1(g), 1(h)], which are leading power effects in 
QCD factorization, are missing in "naive factorization" and "generalized factorization" . 

Readers are referred to Refs. pT^JTSl] for more details. 



C. Chirally enhanced corrections in QCD factorization 



It is observed that QCD factorization is demonstrated only in the strict heavy quark 
limit. This means that any generalization of QCD factorization to include or partly 
include power corrections of ^qcD/^b should redemonstrate the factorization. There are 
a variety of sources which may contribute to power corrections in l/m?,; examples are 
higher twist distribution amplitudes, transverse momenta of quarks in the light meson, 
annihilation diagrams, etc. Unfortunately, there is no known systematic way to evaluate 
these power corrections in general for exclusive decays. Moreover, factorization might 
break down when these power corrections, for instance transverse momenta effects, are 
considered. This indicates that one might have to give up the ambitious plan where all 
power corrections are, at least in principle, incorporated into QCD factorization order by 
order. One might argue that power corrections in B decays are numerically unimportant 
because these corrections are proportional to a small number KqcD/^h — 1/15. But it 
is not true. For instance, the contributions of operator Qq to decay amplitudes would 
formally vanish in the strict heavy quark limit. However it is numerically very important 
in penguin-dominant B rare decays, such as the interesting channels B irK, etc. This 
is because Qq is always multiplied by a formally power suppressed but chirally enhanced 
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factor = 2mp/(m;,(mi + ^2)) ~ C>{1), where mi and m2 are current quark masses. So 
power suppression might probably fail at least in this case. Therefore phenomenological 
applicability of QCD factorization in B rare decays requires at least a consistent inclusion 
of chirally enhanced corrections. 

The chirally enhanced corrections arise from twist-3 light-cone distribution amplitudes, 
generally called 4>p{x) and (pa{x). For light pseudoscalar mesons, they are defined as 



(P(p')|g(i/)^75g(x)|0) = fpfxp [ du e^(^^'-^+^p'-^'Vp(«), (6) 

Jo 

{P{pMyWul,qix)\0) = tfpMP.Zu - p'.z,) du , (7) 

where /ip = ^^^^^ , z = y — x, and mi and m2 are the corresponding current quark masses. 
So, when chirally enhanced corrections are concerned, the final light mesons should be 
described by leading twist and twist-3 distribution amplitudes [p9| : 

u) - /ip75 0p(u) - a^yp'^z- 



{P{p')\Uy)<is{xm , , 

4 JO 



6 



Sa 



Thus it is crucial to show that factorization really holds when considering twist-3 distri- 
bution amplitudes. The most difficult part is to demonstrate the infrared finiteness of the 
hard scattering kernels T/. For more technical details of this proof, readers are referred 
to the literature [ll5| , p!6 . 



With the effective Hamiltonian and QCD factorization formula, B — > P1P2 decay 
amplitudes in QCD factorization can be written as 

Gi. (9) 



AiB P,P,) = ^ E E vX.{PiP2mB) 

V ^ p=u,c j=l,lo 



where Vp is the CKM factor, {PiP2\Qi\Q)F is the factorized matrix element and is the 
same as the definition of the BSW Largrangian ||^. The explicit expressions for the decay 
amplitudes of P — > P1P2 are listed in the Appendixes of |[TT|JT^ . 

For illustration, we give the explicit expressions of af (i = 1 to 10) for B ^ (using 
symmetric LCDAs of the pion). It is easy to generalize these formulas to the case that 
the final states are other light pseudoscalars. Furthermore, we take only some of QED 
corrections into account in our final formula, in particular the QED penguin insertions. 
Now af for P — >■ tttt in NDR 75 scheme is listed as follows: 

+ § + (11) 
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+ E (C4 + Ce + ^eqCg + ^eqCio) log + G{s,)) + GsCsg] , (13) 

g=uAs,c,b 2 2 6 nib J 

+ E (C3 - y)((1 + ^^,)log ^-\-\a<, + G'{8,) + G'(s,)) 

q=d,b " 

+ E (C^4 + Ce + ^eqCs + ^eqCio) ((1 + ^^.) log + + G<^(s,)) 

g=n,d,s,c,fe 2 2 V 3 mfe y 

+(^ + A^)Csg} , (15) 
a. = C. + | + ^fc3(-F-12), (16) 



N An N 

{{C2 + §)((! + log ^ - ^ - + G'(.,) + C^is,)) 



+(C4 + §) E ^e,((l + Ia^) log - ^ - + G'{s,) + C^is,)) 

q=d,b ^ 



+iC3 + j^ + C, + -^) E J {1 + - A^) log ^ + G'{s,) + G'^ is,)) 

q=u,d,s,c,b ^ b 

+(1 + Ia^)C7,] , (17) 
a.^C. + ^ + ^fc,oF, (18) 
«?o = ^^o.|.S^C..-^{(C..|)(llogA,^(.^)_^) 

+(^4 + §)E^e4logA + G(.,)-^) 

+(^3 + § + C^5 + §) E l%il^og^ + G{s,)) + ^GsCr,]. (19) 

q=u,d,s,c,b " 

Here = 3 is the color number, Ci? = (A^^ — 1)/2A^ is the color factor, Sq = m^/ml, and 
we define the other symbols in the above expressions as 

F = -121n-^ - 18 + / + /^ (20) 
rub 

1 1 
f^l dx g{x)<j>{x), Gs^l dxGs{x)(j>{x), (21) 


1 

G{s) ^ J dx G{s,x)(f){x), (22) 
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G'{s) = J dx G'{s,x)(j)p{x), 



1 

(f)^{x) 



dx G (s,x)-—— r, Af, 
6(1 — X) 



dx 



Q(l-x) 



(23) 
(24) 



where 0(x) [0p(x), (pa{x)] is leading twist (twist-3) LCDA of the ejected pion, and the 
hard-scattering functions are 



1 — 2x 

g(x) = 3 Inx — SiiT, Gs{x) 

1 — X 



1 — X 



G{s, x) = —4: J du u{l — u) ln(s — u{l — u){l — x) — ie), 



1 

G'(s, x) = —3 J du u{l — u) ln(s — ^(1 — u){l — x) — ie), 



1 

G"{s, x) = —2 j du u{l — u) ln(s — u{l — u){l — x) — ie) 



u''{l-uf{l-x) 



(25) 
(26) 

(27) 



du 



s — u{l — u){l — x) — ie 



(28) 



The contributions from the hard spectator scattering [Figs. 1(g), 1(h)] are reduced to the 
factor /^^: 



/ 



N Ff-*^(0)m| 



, , di^l dx'^^l dy 



X 



<f>{y) ^ My) 



l-y MB6(l-y)2^ 



• (29) 



There is a divergent integral in /^^. When we do numerical calculation in this work. 



we will simply parameterize this divergence as what was done by Beneke et al. |14 



Xh= dy/y = \n{mi,/ Aqcd) + Que 



l<pH 



(30) 



where (pn is an arbitrary phase, 0° < 360°, and qh is varied from to 3 (realistic) 

or 6 (conservative). Q 



^This parameterization is a little bit different from that in the recent paper ||T^ by Beneke et 
al, but the variation of these two parameterization are almost the same. In addition, we pointed 
out in 1 15 1 that /^^ might be convergent if the transverse momentum kx of the parton is taken 
into account and the Sudakov resummation is invoked. In Ref. |30|, the authors find that value 
of /^^ is really in the range of the above parameterization when the Sudakov suppression is taken 
into account. In recent literature of the perturbative QCD (PQCD) method, the authors 
point out that the threshold resummation can also help to suppress the endpoint singularities, 
especially for the twist-3 level. However, this treatment contradicts the claim of Descotes-Genon 
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We illustrate numerically the scale dependence of af in Table ||. Here we use the 
asymptotic form of the LCDAs of the light pseudoscalar meson, which are 

= 6x(l — x), (31) 
= 1, (32) 
4>ct{x) = 6x{l — x), (33) 

and we set /^^ = in computation because of its uncertainty. 

From the expressions and numerical results of the coefficients a^, some general obser- 
vations are listed below. 

(1) tti (i = 1 — 5, 7, 9, 10) and a6,8'^x invariant so they guarantee a full decay 

amplitude independent of the renormalization scale. We have shown that the coefficients 
ttj (z = 1 — 5, 7, 9, 10) are renormalization scale independent at the order of with 
asymptotic LCDAs in | |15|| , 

"^"'^ (z^6,8), 



dln/x 



However, and ag must be scale dependent because the factorized matrix elements 
multiplied by and Og are scale-dependent. But, 

^ = (. = 6.8), 
dm/i 

This point can also be seen roughly from numerical results in Table |I[ It should be noted 
that the imaginary part of is at the order of a^, and its remaining scale dependence 
must be cancelled by the radiative corrections from higher order of a^. 

(2) The coefficients are gauge invariant. At the leading power of AQcn/nT^b, the 
light meson is represented as its leading twist light-cone wave function (LCWF). In the 
calculation, the Dirac structure of the leading twist LCWF guarantees that the light meson 
can be taken as a pair of coUinear on-shell massless quark and antiquark. Therefore, the 
on-shell condition guarantees the gauge invariance of (z ^ 6, 8). In Ref. [0], the authors 
have shown that at the twist-3 level the light meson can be also taken as a pair of on-shell 
massless quark and anti-quark when the asymptotic form of twist-3 LCDAs are taken. 
Then, a6,8 are also independent of the gauge choice. 

(3) The quantities of ai, a2 receive a special attention, because they are related to tree 
operators Qi,2, so they are numerically large. QCD penguins effects are incorporated to 
the coefficients as-e, in which 03^5 are smaller than 04^6; and 03 -|- 05 is C(10~^). 04^6 have 
the same order as 02, so they will play a dominant role when the effects of tree operators are 
CKM-suppressed, for example in the 6— >s transition. The electroweak penguin coefficients 
ay-Oio are even smaller, og is the largest one among them, and numerically comparable 
with the QCD penguin coefficients 03^5. Thus, ag may play an important role in some 
decay modes, such as 5 ^ n'^K. 



and Sachrajda. They conclude in Ref. |26| that it is impossible to make reliable predictions for 
power corrections to the amplitudes for exclusive two-body B decays which have end-point 
singularities, even if the Sudakove effects are taken into account. We think that this point of 
view needs further careful investigation. 
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(4) The imaginary part of Oj arises from the radiative corrections. Unhke in the gener- 
ahzed factorization, the imaginary part of a, comes not only from the penguin corrections 
(BSS mechanism), but also from the vertex corrections. However, it is suppressed by 

and generally numerically small. Thus, the strong phases evaluated from the BBNS 
approach are generally small. But in some special cases the situation will be different. 
For example, the imaginary parts of 02,4,65 especially for 02, are larger than those of the 
other coefficients, so, when their contributions are dominant in some decay modes, strong 
interaction phases might be large, which would result in a large CP violation. 



D. Contributions of annihilation amplitudes 



Annihilation contributions appear in almost all charmless decay modes B —>■ PP. In 
some cases, they may be important in spite of the power suppression. As emphasized in 



the PQCD method p5|j34 , weak annihilation can give a large imaginary part to the decay 
amplitudes. Then, within the framework of the PQCD method, large CP violations are 
expected. In this work, we will follow Beneke et al. [0, writing the annihilation contri- 
butions in terms of convolutions of "hard-scattering" kernels with LCDAs by ignoring the 
soft endpoint divergences. This kind of treatment is obviously not self-consistent; how- 
ever, it can help us to give an estimation of the annihilation. The corresponding diagrams 
of weak annihilation are depicted in Fig. 2. Here, we will follow the convention in Ref. 



A''''''{B^P^P2)^fBfpJp,T.^^ r h (34) 

The decay amplitudes of weak annihilation for B —>■ PP are listed in the Appendix. 
Similar to the treatment for the endpoint divergence in hard spectator scattering, we 
also parameterize the endpoint divergence in the weak annihilation as follows: 

Xa = I dy/y = ln(m,/AQCD) + 0Ae"^\ (35) 

where (pA is an arbitrary phase, 0° < 0^ < 360°, and is varied from to 3 (realistic) 
or 6 (conservative). It should be noted that there is no correlation between Xa and X^ 
for their different origins. 

We will apply the above formulae to revisit rare B decays in what follows. 



III. INPUT PARAMETERS 



The decay amplitude for B^PP can be expressed by various parameters, such as the 
CKM matrix elements, form factors, Wilson coefficients Cj(/i), LCDAs, and so on. The 
values of these parameters will affect our predictions for CP-averaged branching ratios 
and CP-violating asymmetries. Now we will specify them for use in calculation. 



A. CKM matrix elements 



The CKM matrix in the Wolfenstein parameterization reads 



12 



/ 1-AV2 A AX%p-ir])\ 

VcKM=\ -A 1-AV2 \+0{\^). (36) 

\AX\1- p-tr]) -AA2 1 ) 

Here we take the Wolfenstein parameters from the fit of Ciuchini et ah ||33|: A = 
0.819±0.040, A = 0.2237±0.0033, p = p{l - \^ /2) = 0.224±0.038, rj = r]{l - = 
0.317±0.040 and 7 = (54.8±6.2)°. Correspondingly, we have p = 0.230±0.039, t] = 
0.325±0.039, and y/p^~Trp = 0.398±0.040. If not stated otherwise, we shall use the 
central values as the default values. 



B. Form factors and decay constants 

Decay constants and heavy-to-light form factors are defined by following current ma- 
trix elements: 

(P(g) 1^17^75^2 1 0) = -z/pg^, (37) 



{Piq)\q^,b\B) 



2 2 



+ 



— m] 



■(p-g)/,Fo((p-g2). 



(38) 



For the decay constants we use /^r = 131 MeV, fx = 160 MeV, and fs = 180 ± 40 
MeV. As to the decay constants related to the t] and 77', we shall take the convention in 
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{0\qri5q\v^'\p)) = ^f'ioP 



1 



and 



(0|s75s|?7^^(p)) = -I — '- ^ '— 



{0\uj,u\v^'\p)) = {0\d^,d\v^'\p)) = -^(0|s753|r/W(p)). 

Here we shall not consider the charm quark content in r]^'\ The relations between 77-77' 
mixing and rjs-rjo [the SU{3) octet and singlet] are 

|?7) = cos 6*81778) -sin6'o|r7o), W) = sin 6*81778) + cos 6'o|?7o), (39) 



cu fs 



(40) 



f" = - 

J n 



2-^ cos 6^8 



A 



sin 6*0, 



-2-^ sin 6^8 + cos 6*0, 



(41) 
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Here /§, /o are decay constants of r^g, rjQ respectively, 6^8, Oq are the mixing angles. Their 
values are /g ~ 1.28/^ ^ 168 MeV, /o ^ 1.20/^ ^ 157 MeV, Os ^ -22.2°, and 
^0 ~ -9.1°. 

We take F^"" = 0.28±0.05, which is obtained by the light -cone sum rule [PH]] . For 
B^K, 7]^''^ transitions, the following relations are applied to give their corresponding 
values: 



C. Quark masses 



There are two types of quark mass in our analysis. One type is the pole mass which ap- 
pears in the loop integration. Accompanied by the 0{aj) corrections to the hadronic ma- 
trix elements, the pole masses contribute to coefficients through the functions G{s,x), 
G'{s,x) and G"{s,x), or G{s), G'{s) and G"^(s). Here we fix them as 

m„ = = rris = 0, = 1.45GeV, rrif, = 4.6GeV. 

The other type quark mass appears in the hadronic matrix elements and the chirally 
enhanced factor through the equations of motion. They are renormalization scale 
dependent. And their values have been summarized in [^]. Here we shall use the 2000 
Particle Data Group (PDG2000) data for discussion: 

mb(rnb) = 4.0 ~ 4.4GeV, m,(2GeV) = 75 ~ 170MeV, 



mrf(2GeV) = 3 ~ 9MeV, m„(2GeV) = 1 ~ 5MeV. 

Here, we take the central values of the fequark and s quark as default values. Because the 
current masses of light quarks are difficult to fix, we would like to take 



fu 



JrjC) 



We think that it is a good approximation in our calculation. Using the renormalization 
group equation, we can get the corresponding running at the scale /i = (9(mb): 

r^{mb/2) = 0.85, r^{mb) = 1.14, r^{2mb) = 1.42. 



IV. NUMERICAL ANALYSIS 
A. Branching ratios and comparison with data 

In the B meson rest frame, we have two-body decay width 

TiB^P.P^) = ^%A{B^P,P2)\\ (43) 
87r 
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where 



bl 



— (mpj + mp^y m\ — {mp^ — mp^Y 



2m B 

The corresponding branching ratio is 

BR{B^P,P,) = £(^Z!^, (44) 

J- tot 

where r = 1/Ttot- In our calculation, we take r = 1.548x10^^^ sec for , and r = 
1.653x10"^^ sec for B^. In addition, because we work in the heavy quark limit, we take 
the mass of light meson as zero in the computation of phase space; then \p\ = tjib/^. 

We show experimental data and theoretical predictions for the CP-averaged branching 
ratios for B^PP in Table |TTT| and Table ^V\, respectively. In Table |V|, we give the averages 
of the branching ratios of B—>-PP and their CP-conjugate modes with the default values 
of the parameters at three different renormalization scales fi = mb/2, rrib, 2mb. From 
Table |IT| and Table |I^ some remarks are in order. 

(1) For the decay modes B — > tttt and B nK, our predictions are consistent with 



those by Beneke et al. |jT6|. In these modes, with the default values of the parameters, we 
find that the CP-averaged branching ratios of B^ tt^ti^ (~ pure tree) and B^ K^n^ 
(pure penguin) are in good agreement with the experimental measurements. Therefore 
these two decay modes can be considered as good probes of QCD factorization for B 



meson decays to two light mesons [^]. For the other modes B — >■ tttt and B — > tcK, 
the predictions of the CP-averaged branching ratios obtained by using the default values 
of the parameters seem not very good compared with the experimental measurements. 
In particular, the theoretical prediction for the branching ratio of 5° tt+tt^ is larger 
than the experimental data. The same situation also occurs in generalized factorization 



rT| , p!2[| . For the other decays modes B — >■ tiK, our predictions seem smaller than the 
measurements. But taking the uncertainties of the input parameters into account, such 
as the CKM matrix elements (especially the angle 7 and iKiftI), the endpoint divergences in 
hard spectator scattering and weak annihilations, form factors, etc., it is possible that the 
theoretical predictions from the BBNS approach can be in accord with the experimental 
results. In Refs. Il37|j36[] , the author gave a global fit to the CP-averaged branching ratios 
of — > TTTT, nK within the framework of the QCD factorization. With the best fit values, 
the branching ratios of B —>■ mr and irK are in very good agreement with the experimental 
measurements very well. [Although, in the best fit, value of 7 is slightly too large (about 
90°) and \Vub\ is comparably small, they are still compatible with the standard fit using 
semileptonic decays, K — K mixing and B — B mixing.] 

(2) We note that experimental data for CP-averaged branching ratios for B^Krj' 
are about four times larger than our theoretical predictions with default values of the 
parameters. The results under the generalized factorization (GF) framework are similar 
|TT| , p!2[| . In Refs. |jll],|l2|, the charm quark content in r]' is considered; however, this mech- 



anism cannot give a good explanation for experimental data yet. The unexpectedly large 
data have triggered considerable theoretical interest in understanding the mechanism of 
77' production, and have been widely discussed in the literature So far, it is gener- 
ally assumed that digluon fusion mechanism could enhance rj' production, but the shape 
of form factor of the vertex g*g*r]' brings a very large uncertainty to predictions of the 
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contribution from digluon mechanism Recently, M.Z. Yang and Y.D. Yang gave a 
calculation for B rj^'^'P within the framework of the BBNS approach [^], in which the 
digluon mechanism is considered. In that paper, the authors compute the vertex of g*g*r]' 
in perturbative QCD, and find the branching ratios of B ^ Kr]' are really enhanced and 
in agreement with the experimental measurement. But, as mentioned in our Introduction, 
we think that the consistency of their perturbative calculation is questionable due to the 
endpoint behavior. 

(3) For the decays B 

± ^ 7r^n^'\ K^T], the CP -averaged branching ratios predicted by 
the BBNS approach are at the order of 10~^. We expect that these decay modes can be 
observed soon at the BaBar and Belle experiments 0. 

(4) For the decays B^ — > vr^vr", '7i^ri^'\ ri^'^ri^'\ our predictions show that their CP- 
averaged branching ratios are very small and about 0(1O~^) ~ O{10~^). If these predic- 
tions are reliable, it is very hard to observe these decay modes at the BaBar and Belle 
experiments. 

(5) In the BBNS approach, the predictions for the branching ratios of -B ^ KK are 
highly suppressed, especially for B^ —>■ K~^K~ , which is purely from the weak annihilation 
contributions. In Ref. ||4^, the authors evaluate the decays B — > KK by use of the PQCD 
method, and also found that their branching ratios are very small. These decay modes 
can be a good probe for long-distance final state interaction (FSI). On the other hand, as 
we show in the Appendix, if the endpoint divergence is universal, the measurement 
for the branching ratio of B^ K^K~ may give a constraint for the range of Xa- 

In Table [rV| , for comparison, we also give the estimations for the CP-averaged branch- 
ing ratios oi B PP using the naive factorization (NF) approximation. Again, we list 
the predictions by the BBNS approach with and without considering the weak annihilation 
also. From these results, some general observations are given as follows. 

Obviously, predictions on branching ratios for B—^PP using the BBNS approach are 
less scale dependent compared with those in the NF framework, especially in the b—>-s 
transition, such as B^Ktt, Kr]^'\ Noting that are calculated at one-loop level, if 
the corrections from high order of as were taken into account, the scale dependence of 
theoretical predictions might be further reduced. 

Generally, the contributions from weak annihilation are highly power suppressed. How- 
ever, due to the endpoint divergences, they may give observable contributions in our re- 
sults. In B ^ ttK, we can see that the weak annihilation can give about 5% ~ 10% 
enhancement. Note that the results of Table ^ are computed by the default parame- 
ters. If we varied the parameters for the endpoint divergences in weak annihilations, it 
would give much larger uncertainties to our predictions. This will be shown in the next 
subsection in detail. 



^Recently, Belle and BaBar collaborations report their measurements of the branching ratio of 
the decay n^-q' which is consistent with zero. The 90% CL upper hmit is BR(5± 

■K^r]') < 7 X 10-6 (at Belle) or BR(S± tt^t]') < 12 x 10"^ (at BaBar). 
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B. Uncertainty of predictions 



In Fig. 3, based on the BBNS approach, we show the dependence of branching ratios 
for B — i> 7171, 7tK and Ki]' on the weak phase 7, scanning all ranges of the variation 
of parameters, including the CKM matrix elements, form factors, endpoint divergence 
J^dx/x, and so on. Note that weak annihilation effects do not contribute to -B^ — >■ tt^tt". 
We can see that our theoretical predictions give very wide dot shades for the branching 
ratios within the range of parameters. So it is possible that the predictions for CP- 
averaged branching ratios of the decays B ^ 7171, 71 K are simultaneously consistent with 
the experimental data by use of the same input parameters. The results of the global fit 
in Refs. [[I6| , |36| , p7| show that, with the best fit values of the input parameters and within 
the framework of QCD factorization, the predictions for the branching ratios of B — vrvr 
and 7iK can fit the experimental data simultaneously with good l^dof- 

In Fig. 4, we show the dependence of the CP-averaged branching ratios of B^ 7[^K^, 
B^ — > K^7T^ and B^ vr+vr" on the CKM matrix elements, form factors, and endpoint 
divergence in hard spectator scattering and weak annihilations. From Fig. 4, some 
remarks are given in order. 

For the decays B^ —>■ 7t^K^ and B^ —>■ K^7f^, the CP-averaged branching ratios 
weakly depend on |\4fe/Vcf,| because that the decays B nK are dominated by penguin 
diagrams and \vu/vc\ = \(y*gVub) / {V*gVcb)\ is highly suppressed by (9(A^) in the b ^ s 
transition. For the case of B^ vr^vr", although it is a tree-dominant (6 u) process, 
the penguin pollution plays an important role. For illustration, the decay amplitude of 
B^ TT+TT^ reads 



a1 



+ a^o + (a^ + 4)rJ}. (45) 



Note that \vu\ = \Vud^ub\ is at the same order as \vc\ = \V*^Vcb\ and \{ai + aQr^)/ai\ ~ 0.1, 
so the branching ratio of 5° — >• tt+tt" is sensitive to both |Ktfe/Kf,| and 7 = arg(V^*j,). 
(In contrast, the decay amplitude of B suffers less pollution from the penguin 

contribution, which is proportional to the small electroweak penguin coefficient ag + aio — 
ay + asr^.) Therefore, the variation of |Vu6/Kf,| gives large uncertainty to the prediction 
for the branching ratio of B^ — »• tt+tt". 

All the three modes are dominated by the form factor F^'^{0). Varying F^'^(O) from 
0.23 to 0.33, it gives about ±35% uncertainty to the branching ratios. By using the ratios 
of the branching ratios, the uncertainties caused by the form factor can be largely reduced. 
This will be shown later. 

The endpoint divergence in the hard spectator scattering does not bring large uncer- 
tainty to our predictions. The reason is that the endpoint divergence contains only a 
single logarithmic term and its coefficient is suppressed by and the color number N. 

Obviously, the endpoint divergence in weak annihilation causes large uncertainty, in 



particular, to the decays B —> nK. In Ref. |lg], the authors also showed this point. The 



reason is that the endpoint divergence in weak annihilation is dominated by a product of 
two logarithmic terms. Within the range of the parameters g and 0, our predictions by 
the BBNS approach are in agreement with the experimental measurements. 

Therefore, we find that weak annihilation is the main source of the uncertainties in 
the BBNS approach. To obtain a more precise prediction, further analysis is needed of 
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the weak annihilations. On the other hand, there are other possible power suppressed 
corrections in the decays B — > PP. Some researchers pointed out that the elastic and 
inelastic long-distance FSI can help us obtain small branching ratio for tt+tt^ at 

small 7 [JT,i2|. However, with the elastic FSI, in order to get satisfying results, one must 
introduce a large strong phase difference which is unbelievable in the heavy quark limit. 
In particular, in Ref. the author takes the inelastic rescattering B — > DD — > vrvr 
into account and finds that the small branching ratio of — vrvr can be obtained at 



small 7 without fine-tuning the input parameters. Furthermore, in Ref. Isola et 
al. give detailed research on the nonperturbative charming penguin effects on i? — > tiK 
and the predictions are in better agreement with the experimental data. However, all 
treatments to these power corrections are model-dependent. Up to now, there is no 
systematic way to deal with these complicated nonperturbative interactions. Maybe the 
precise measurement of the CP asymmetries for the decays B PP is a chance to 
distinguish these models. 

Here we emphasize that the decay amplitude of 5° K~K^ comes completely from 
the effects of weak annihilations under the BBNS approach, because contributions of weak 
annihilations are supposed to be power suppressed, so the CP-averaged branching ratio for 
5 (5°) ^ K+K- is small in Fig. 3. At the same time, we can also see that uncertainties 
of theoretical predictions for these decay modes are very large (the branching ratios range 
from 10~^ to 10"''). As mentioned above, we expect that we can extract some useful 
information about the endpoint divergence Xa from B^ K^K~ . However, the present 
experimental upper limit is much greater than our predictions within the range of the 
parameters. If this decay mode can be measured precisely in the future, then we will be 
able to obtain some useful information on weak annihilations and FSI. Consequently, the 
endpoint divergence Xa can be determined effectively. 



C. Ratios of the branching ratios 

As mentioned above, the decay rates for B^PP depend on a number of parameters 
that cannot be predicted with the BBNS approach itself. Thus, these parameters result 
in ambiguity of the predictions for the CP-averaged branching ratios. In order to reduce 
uncertainties from the renormalization scale and weak annihilation effects, etc., as much 
as possible, the ratios of CP-averaged branching ratios are good quantities. Furthermore, 
many methods have been suggested to give bounds on 7 using decay modes of B^nn, nK, 
in which the ratios of branching ratios play important roles. We hope that our calculation 
can give some constraints on 7. In Fig. 5, we show some ratios of CP-averaged branching 
ratios versus the weak phase 7 (note that we use PDG2000 data for tb+/tbo = 1.062), 
where the slashed-line bands correspond to experimental data listed in Table ^ with one 
standard deviation (la). 

Obviously, using the ratios of branching ratios, and comparing with Fig. 3, the width 
of dot-shades becomes narrow due to the reduction of the uncertainty. From Fig. 5(b), the 
ratio BR(S*^ K^) / {2'BIi{B^ n^K^)) can be made consistent with the experimen- 
tal data at small 7. However, the ratio rB+Zr^o [BR(S° 7r+7r")/(2BR(5° tc^tt^))] 
tends to be consistent with the experimental data at large 7 in Fig. 5(e). Considering 
the correlation between these quantities, a scientific choice is to give a global fit for these 



branching ratios and/or the ratios of the branching ratios. In Ref. [|T6l, the authors give a 
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global fit for (p, rj) from B rare hadronic decays with the QCD factorization. The result 
favors large 7 (about 90°) and/or small \Vub\ 



There are many discussions on ratios of branching ratios, for instance in [0,0) in 
order to determine coefficients a^, and extract the CKM matrix parameters, etc. In 
particular, the decay modes B^Kir are widely discussed to get information on 7; readers 
interested in these topics are referred to, for instance, pi| , ^ 
our theoretical results show that 



From Figs. 5(a) and 5(b), 



2BR(5±-^7rOir±) _ BR(50^7rTi^±) 
BR(5±^7r±K0) ~2BR(50^7rOiro) 



(46) 



This approximate relation can also be obtained from the isospin decomposition of the 
effective Hamiltonian ^ |^6|. The experimental data listed in Table |T| give 



2BR(7rOir±) 
BR(7r±fs:o) 



1.18±0.34(BaBar), 



2BR(7rOi^o) 



1.02±0.41(BaBar), (47) 



which are consistent with Eq.(E6| 



D. CP- violating asymmetries 

CP violation in B meson decays can be either direct or indirect. Typically, direct CP- 
violating asymmetry arises when the tree and penguin amplitudes interfere. CP violations 
in decays of charged B meson are purely direct; as these decays are all self-tagging, they 
are the definitive signals of CP violation if seen. The direct CP-violating symmetry is 
defined by 

-^^^ - V{B-^f-) + V{B+^n- ^^^^ 

Because the fiavor eigenstates 5^ and B^ mix due to weak interactions, direct and in- 
direct CP violation occur simultaneously in neutral B meson decays, and time-dependent 
measurements of CP-violating asymmetries are needed. 

. r(gW-/)-r(B°W-/). («) 
r(B (i)^/) + r(B»W^/) 

For B^^f (but B^y^f) and B^-^f (but B^-/^f), the CP- violating asymmetries are similar 
to those of charged B mesons. For B -^(/ = /, CP\f) = ±\f))^B^ , i.e., the final states 
are CP eigenstates, the time-integrated CP-violating asymmetries are: 



^In Ref. H, tiie isospin amplitudes are defined as B1/2 = ^273(1/2, ±1/2|Wa/=o|1/2, ±1/2), 
which arises from QCD penguin operators, while vli/2 = ±-\/2/3(l/2, ±1/2|'Ha/=i|1/2, ±1/2), 
^3/2 = ^1/3(3/2, ±1/2|'Ha/=i|1/2, ±1/2), which include contributions from tree and 
electroweak operators. Moreover, \Bi/2\^0{W-^), 1^1/21^^0(10-^), and \As/2\^0{W-^). 
The decay amplitudes can be described as A{B^^7r~^K^) = VI3/2 ± ^1/2 + -^1/2) ^i^d 
\/2A{B^^Tr'^K^) = 2A^/2 — ^1/2 — B1/2. Neglecting effects of ^3/2 and ^1/2, we can get 
a good approximation to Eq.(|4^. 
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1 X 

+ \Xfr " 1 + iA,p ' - v,,v:, ^ ^ 

with the PDG2000 data value x = Am/T = 0.73±0.03 for B°-B° system. Here a,/ and 
a^+e' are due to direct and mixing-induced CP violation, respectively. 

Recent experimental measurements of CP-violating asymmetries Acp for B^ttK, 
Ki]' and a^', a^+e' for n^n^ are listed in Table 0. These results contain very large 

errors and they are consistent with zero. Precise measurements of direct CP violation 
and time-dependent CP violation are expected in the near future. 

Theoretically, CP-violating asymmetries will not be very large in the BBNS approach 
in principle, because the strong phases are suppressed by as or AQcn/^b- However, from 
Table || we know that 02,4,6 have large imaginary parts. When they are dominant in some 
decay modes, large CP violations are expected, for instance, in B (i?°)^7r°7r*', TT^ri^'\ 

With the BBNS approach, CP-violating asymmetries depend on several variables, such 
as the form factors, CKM parameters, the scale /x, and so on. Our investigation indicates 
that CP-violating asymmetries for all B^PP depend very weakly on form factors; the 
same conclusion can be obtained with the GF approach. So we show results only with 
the default values of the form factors in the following calculations. As for the dependence 
of CP-violating asymmetries on scale /i, there are six decay modes B (i?°)^7r°7r°, 7i^t]^'\ 
j^iOrqi') singled out because they are strongly fi dependent compared with other decay 
modes. In these six decay modes, 02 is involved and not CKM-suppressed, and it is sensi- 
tive to the renormalization scale. This point has been shown in Table |I|. In particular, its 
imaginary part is larger than the real part, which might result in large CP violations in 
these decay modes as we show in Table ^I[ We also list CP-violating asymmetries versus 
the weak phase 7 in Table |Vli| and Table |Vlli]. The endpoint divergence can cause large 



theoretical uncertainties; its influence on CP- violating asymmetries is shown in Table IX. 



From Table |VT| to Table IX, we can clearly see that theoretical predictions on CP- 



violating asymmetries are compatible with measurements within one standard deviation 
(Icr). Unfortunately, the present experimental uncertainties are too large to draw any 
meaningful conclusion. In addition, the weak annihilations have great influence on CP- 
violating asymmetries, so their contributions to decay amplitudes must be included in 
estimation of CP violation in B meson decays with the BBNS approach. Moreover, the 
power corrections in order of KqcD/^b which are not included in the QCDF master 
formula (^). This needs further investigation. 

Here we would like to point out that 5— ^vrvr are usually considered to determine the 



CKM angles a and 7 using CP-violating asymmetry [^,@ . From Table |VII| and Table 
Villi , we can see that Aiji^n^), especially the mixing-induced a^^^iiVo) term, are sensitive 
to the CKM angle 7. This mode is strongly polluted by the penguins. As to the large 
penguin pollution, recently, Fleischer suggested a new approach Using U-spin SU(3) 



flavour symmetry, 5^— ^-vr+vr" and Bg^K^ K~ are related to each other by interchanging 
the d quark and s quark, which possibly allows a simultaneous determination j3 and 7. 
His analysis gives 7 = 76° and 0^ = 2/5 = 53° without considering the U-spin-breaking 
effects. 
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As we see from Table [VI I| and Table [VIII] , time-integrated CP- violating asymmetries 
Acp for B ^KgTc, K^r]^'^ are completely dominated by mixing-induced a^+e' term. But 
we still do not know the i]^'^ production mechanism. Hence it is useless to discuss the CP- 
violating asymmetries Acp{Kgr]^'^). As to CP-violating asymmetries for B ^K'^K~^ 
they do not contain the direct CP-violating a'^ term, and arise only from 0^+^/ term. 
Moreover, they are sensitive to the CKM angle 7. However, we know that with the BBNS 
approach the branching ratios for neutral B meson decays into K'^K~ have only weak 
annihilations contributions, which are scale 11 dependent and contain large theoretical 
uncertainties due to their soft nonperturbative effects. 



V. CONCLUSIONS 



(1) In the heavy quark limit, neglecting the effects of order Agc£,/m6, the contributions 
that are non-factorizable under the OF framework can be factorized with the BBNS 
approach, and they are perturbatively calculable from first principles, at least at order of 
as- The BBNS approach provides decay amplitudes with strong phases which are usually 
small, at order of as and/or KQcn/T^b- 

(2) CP-averaged branching ratios of B-^PP in the BBNS approach are less renormal- 
ization scale dependent than those under the NF framework, especially when the chirally 
enhanced corrections are taken into account; generally, with the appropriate parameters, 
predictions with the BBNS approach are consistent with present experimental data. How- 
ever, the theoretical predictions include very large uncertainties which arise from various 
parameters, such as the form factors, CKM matrix elements and soft endpoint divergence, 
etc. Our results show that the theoretical uncertainties from the endpoint divergence Xa 
are considerable. Likewise, CP-violating asymmetries for B^PP in the BBNS approach 
are also sensitive to variation of endpoint divergences. In particular, decay amplitudes for 
B {B°)^K+K- arise completely from weak annihilations with the BBNS approach; we 
can obtain very useful information on the effects of FSI and soft weak annihilations with 
precise experimental data in the future. 

(3) The present experimental data for non-leptonic B meson decays are not precise 
enough yet, especially for the measurements of CP violation. On the other hand, there 
are still many uncertainties in the theoretical framework, for instance, the hard spectator 
scattering, weak annihilations and other potential power corrections. Great advances in 
both experiment and theory in the near future are expected to determine the CKM angles 
precisely. 
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APPENDIX: THE ANNIHILATION AMPLITUDES FOR B^PP 
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(A2) 
(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 



(AlO) 



(All) 
(A12) 

(A13) 

(A14) 
(A15) 
(A16) 



The annihilation coefficients (61, 62), (^3, ^4) and (63"', If^) correspond to the contributions 
from tree, QCD penguins and electroweak penguins operators, respectively. They are 
related to final state mesons. Using the asymptotic light cone distribution amplitudes of 
the mesons, and assuming 5'?7(3) flavor symmetry, they can be expressed as [F* 

^^\c^A' + Af(G^ + N,C^) 



b. = %C,A\ 

c 

b. = %C2A\ 



b, = ^A\c, + C,), 



"a 



c 



Cs^A' + Af (Cj + N^d 



c 



(A17) 
(A18) 
(A19) 
(A20) 
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and 



18 Xa - 4 + — + 2r2X 



-) 

3 ^ 

^12nasr^(2Xl - Xa), 



(A21) 
(A22) 



where Xa = j\dy /y is a logarithmically divergent integral, and will be phenomenologi- 
cally parameterized in the calculation. 
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TABLES 



TABLE L Wilson coefficients in NDR sclieme. Tlie input parameters in numerical calcu- 
lation are fixed: as{mz) = 0.1185, aem{mw) = 1/128, mw = 80.42GeV, mz = 91.188GeV, 
mt = 168.2GeV, = 4.6GeV. 





u = mh/2 

r (J/ 


(J 


H = 2m(, 




1.136 


1.080 


1.044 


Co 


-0.283 


-0.181 


-0.105 




U.UZ 1 


n m A 


n nno 




-0.050 


-0.035 


-0.024 




0.010 


0.009 


0.007 


Ce 


-0.063 


-0.041 


-0.026 


C7/ aem 


-0.020 


-0.004 


0.019 




0.082 


0.052 


0.033 




-1.339 


-1.263 


-1.201 


Cio/Oem 


0.369 


0.253 


0.168 


C77 


-0.341 


-0.304 


-0.272 




-0.160 


-0.145 


-0.132 



TABLE II. 


Numerical values of coefficients Oj 


with = 4.6GeV, rric 


= 1.45GeV and /^^ = 0. 




;[i = m{,/2 


li = mb 


H = 2mh 




1.070+0.0261 


1.046+0.0131 


1.027+0.0061 


02 


-0.019-0.106i 


0.023-0.0791 


0.062-0.0641 


03(10-^) 


93.662+46.5851 


73.837+25.7291 


51.675+14.6461 


<(io-4) 


-330.78-179.231 


-297.17 -146.841 


-267.83 -125.461 




-403.19-55.9031 


-351.48-54.3371 


-311.41 -51.2261 


a5(10-4) 


-98.679-58.6971 


-67.453-30.1391 


-41.612 -15.867i 


<(io-^) 


-556.39-159.841 


-411.61-136.651 


-322.84 -120.021 


ag(10-4) 


-597.2-34.6521 


-442.22-42.7421 


-347.41 -44.6631 


o^x(10"') 


-474.17-136.221 


-469.29-155.8 1 


-459.82-170.951 


air^(10-4) 


-508.95-29.5311 


-504.19-48.7311 


-494.81-63.6141 


07(10-^) 


0.458+0.5971 


1.248+0.2991 


2.591+0.1571 


a|(10-4) 


6.980-0.661 


4.255-1.1171 


2.336-1.4651 


4(10-^) 


6.884-0.3651 


4.076-0.5661 


2.092-0.7171 


«^x(10"') 


5.949-0.5621 


4.852-1.2741 


3.328-2.0871 


a§r^(10-4) 


5.867-0.3111 


4.647-0.6451 


2.980-1.0211 


09(10-^) 


-97.902-2.6861 


-94.935-1.4531 


-91.732-0.8011 




7.414+9.0871 


3.136+6.1361 


-0.956+4.2611 


afo(10-^) 


7.244+9.3771 


2.817+6.6791 


-1.389+4.9981 
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TABLE III. Experimental data of CP-averaged branching ratios for B^tttt, ttK, Kt]' in unit 
of 10~6. 



Decay Modes 



CLEO 1 1 



BaBar || 









4.1±1.0±0.7 






5.6t|^±1.7 


5.ij:?:^±o.8 






17.2l|^±1.2 


16.7±1.6±1.3 






1 1 (5+3.0+1.4 
11.0_2 7_i 3 


io.8li:9±i.o 






18.2l^;|;±1.6 


18.2l|^±1.7 






14 R+5.9+2.4 
-■-^•0-5.1-3.3 


8.2l|^±l.l 


B-- 




80l^°±7 


70±8±5 


b'- 




89lJ^±9 





Belle g] 



7.8^ 



5.6l2;|i±0.4 

j+3.8+0.8 
-3.2-1.2 
iq 0+3.4+1.5 
iy.o_3 2_o,6 

-iq 7+5.7+1.9 
-4.8-1.8 

-"-O-U.s 9_2.7 

79l}?±9 



TABLE IV. Numerical predictions for CP-averaged branching ratios (in unit of 10 ^) for 
B^PP, in the framework of NF and QCDF, where BR-'' and BR-^+'' denote the CP -averaged 
branching ratios without and with the contributions from weak annihilation, respectively. The 



experimental data is from PDG2000; * denotes uncorrelated averages of Table III 









IJ, = mb/2 




IJ, = mi, 






/i = 2mh 








Decay 


NF 


QCDF 


NF 


QCDF 


NF 


QCDF 


Exp. 


Modes 


BR 


BR-' 


BR^+" 


BR 


BR-' 


BR-'+'' 


BR 


BR^ 


BR-'+'' 




b' 




Q 7/1 


10.06 


10.49 


Q no 
y.uy 


9.67 


9.97 


O.Do 


9.36 


9.58 


A /i_i_n Q* 


B' 




5.58 


5.26 




6.20 


5.36 




6.76 


5.52 




5.7±1.5* 


B- 




3.41 


3.27 


3.22 


3.59 


3.28 


3.26 


3.74 


3.34 


3.33 


< 15 


B- 




2.26 


2.23 


2.18 


2.53 


2.26 


2.24 


2.71 


2.35 


2.34 


<12 i 






0.17 


0.17 


0.18 


0.12 


0.15 


0.16 


0.07 


0.14 


0.15 


< 8 


b' 




0.038 


0.048 


0.054 


0.037 


0.045 


0.053 


0.026 


0.049 


0.057 


< 11 






0.12 


0.11 


0.18 


0.15 


0.099 


0.13 


0.23 


0.10 


0.12 


< 9.3 






0.077 


0.078 


0.104 


0.095 


0.076 


0.093 


0.12 


0.081 


0.093 


< 18 


b' 




0.009 


0.016 


0.022 


0.029 


0.015 


0.022 


0.048 


0.021 


0.027 


< 47 


b' 




0.059 


0.073 


0.066 


0.11 


0.073 


0.073 


0.15 


0.086 


0.09 


< 27 


b" 




10.66 


9.56 


10.45 


6.40 


8.59 


9.13 


3.4 


7.78 


8.09 


17.4±1.5* 


b' 




4.99 


4.52 


4.96 


2.90 


4.01 


4.27 


1.44 


3.57 


3.72 


10.3±2.6* 


B- 




7.65 


7.04 


7.56 


4.83 


6.36 


6.67 


2.82 


5.78 


5.96 


12.1±1.6* 


B- 


- 


14.43 


13.45 


14.48 


8.8 


11.99 


12.61 


4.79 


10.72 


11.08 


17.3±2.5* 


B' 




2.91 


2.67 


2.78 


1.68 


2.31 


2.37 


0.94 


1.96 


2.0 


< 14 


B- 




16.6 


16.72 


17.82 


12.47 


15.74 


16.87 


7.49 


15.75 


16.55 


74.9±6.6* 


b' 


-^K T] 


1.93 


1.70 


1.79 


1.0 


1.45 


1.50 


0.45 


1.21 


1.24 


< 33 




-^A T] 


17.07 


17.33 


18.13 


12.82 


16.26 


17.12 


7.79 


16.14 


16.73 


59.6±9.8* 


B- 




0.78 


0.75 


0.79 


0.47 


0.67 


0.69 


0.26 


0.59 


0.61 


<2.5 i 


b' 




0.73 


0.71 


0.84 


0.44 


0.62 


0.70 


0.24 


0.56 


0.60 


< 17 


b' 


^K+K- 






0.044 






0.033 






0.03 


<2.5 i 
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TABLE V. Experimental data of CP-violating asymmetries Acp for B^ttK, K-q' , 
Ge+e' for decay ^tt+tt" in unit of percent. 



and 



CLEO 



BaBar 



Belle 

A ^+18.6+1.8 
^•^-16.7-2.1 
1+22.2+5.5 
-19.6-1.7 
Q Q+43.0+2.0 
^•"-34.3-6.3 

6± 15± 1 



Acpin'^K^) 
Acpin^K^) 
Acp{K^i) 



-4±16 
-29±23 
18±24 
3±12 



-7±8±2 
0±18±4 
-21±18±3 



-25^1^±14 



-5.9^ 



a£'(7r"'"7r ) 
ae+e'(7r+7r") 



TABLE VL The dependence of time-integrated CP-violating asymmetries a^/, Oe+e', and 



Acp (in units of percent) for B {B^)—>-Tr^Tr^ , Tr^rj^'^ ri^'^ rj^'^ on renormalization scale n, with the 
default values of various parameters in the BBNS approach. 



Decay Modes 


scale iJ, 




f+a 

ay 


"e+e' 


f+a. 


A-^ 






mb/2 


-57.3 


-69.3 


-32.0 


26.5 


-52.6 


-32.6 




nib 


-45.1 


-79.0 


-66.8 


-12.4 


-61.3 


-57.4 




2mb 


-38.3 


-82.4 


-89.2 


-46.3 


-67.5 


-75.8 




mb/2 


35.5 


62.7 


9.72 


-30.7 


27.8 


26.3 


B\B^)^7r% 


rub 


3L5 


63.3 


12.2 


-19.3 


26.4 


32.1 




2m.b 


29.5 


65.6 


15.8 


-9.75 


26.8 


38.1 




mb/2 


44.2 


84.8 


25.9 


-35.3 


41.2 


38.5 


B\BO)^Tr°rj' 


nib 


41.8 


81.0 


19.6 


-28.2 


36.6 


39.4 




2nib 


37.1 


79.3 


18.3 


-17.4 


33.0 


43.5 


B°{B°)^riri 


nib/ 2 


55.5 


47.3 


58.8 


51.4 


64.2 


55.4 


rrib 


43.1 


36.4 


71.1 


69.9 


61.9 


57.0 




2m.b 


34.7 


26.6 


82.6 


84.4 


62.0 


57.6 


B"{B^)^r]'r)' 


mb/2 


30.7 


76.8 


95.2 


52.7 


65.4 


75.2 


nib 


37.8 


65.6 


92.3 


71.9 


68.6 


77.1 




2nib 


36.7 


48.1 


92.9 


87.0 


68.2 


72.8 


-^{B^)^r]r]' 


mb/2 


55.6 


49.5 


74.1 


77.8 


71.5 


69.3 


nib 


46.0 


35.3 


81.7 


87.3 


68.9 


64.6 




2mb 


37.4 


23.5 


88.9 


94.2 


66.8 


60.2 
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TABLE VII. CP-violating asymmetries ae'{%) and ae_|_e/(%) for neutral B meson decays with 
default values of various parameters using the BBNS approach. 



Decay 




7 = 


60° 






7 = 


90° 






7 = 


120° 




Modes 


X 




T 


T-\-a 


1 

4 


+■ 1 — 

a{, 


1 


1 „ 

T-\-a 
4+e' 


7 

4 


t H-ft 

a^, 


T 


T-\-a 




4.17 


11.5 


55.1 


54.8 


5.85 


16.1 


-30.7 


-30.6 


6.64 


17.7 


-92.0 


-91.1 


TT^r] 


31.2 


60.8 


11.4 


-21.1 


25.0 


43.5 


7.27 


-22.1 


16.6 


27.3 


4.21 


-16.0 




41.6 


77.4 


18.4 


-31.8 


34.7 


54.5 


11.5 


-35.8 


23.5 


34.0 


6.54 


-26.3 


ttOttO 


-41.5 


-76.6 


-65.3 


-16.0 


-25.6 


-56.9 


-49.8 


-23.7 


-15.1 


-36.4 


-32.1 


-18.8 


rjr] 


45.7 


38.5 


70.6 


69.9 


53.6 


44.4 


60.2 


61.4 


47.2 


38.4 


43.3 


45.0 


rj'rj' 


/1 1 7 


lU.o 


CO 7 


DO. 4 


DO. ( 


oD.y 


DO.o 


AO O 


oy.y 


Q1 1 

oi.i 


QQ A 


Zl.U 


ijrj' 


49.7 


38.4 


81.0 


87.4 


65.5 


53.4 


68.8 


80.5 


66.1 


58.0 


50.7 


65.7 




1.78 


2.37 


77.0 


76.6 


2.14 


2.83 


74.1 


73.7 


1.92 


2.54 


58.8 


58.3 




3.08 


3.62 


77.8 


77.4 


3.72 


4.36 


75.1 


74.7 


3.38 


3.95 


59.9 


59.3 




-1.92 


-2.20 


73.9 


74.3 


-2.25 


-2.58 


70.3 


70.8 


-1.97 


-2.26 


54.7 


55.2 




21.0 


20.8 


9.08 


3.74 


16.6 


16.1 


6.30 


2.16 


10.9 


10.5 


3.86 


1.17 


K+K- 









72.9 









3.47 









-68.0 



TABLE VIII. CP-violating asymmetries Acp{%) for decays of B—>PP with default values 
of various parameters using the BBNS approach. 



Decay 
modes 


7 = 


60° 


7 = 


90° 


7 = 


120° 








Aj+a 
•^CP 


•^CP 


Al+a 
•^CP 


S'(50)^7r+7r- 


28.9 


33.6 


-10.8 


-4.08 


-39.6 


-31.9 




0.07 




0.09 




0.08 




B^^-K^r] 


11.5 


25.5 


18.6 


40.3 


27.2 


55.4 


B^^-K^T]' 


5.08 


19.3 


7.55 


28.1 


9.16 


32.9 




25.8 


29.7 


19.8 


17.9 


12.8 


10.2 




35.9 


35.4 


28.1 


18.5 


18.5 


9.60 




58.1 


57.6 


40.4 


48.4 


25.2 


32.7 


B°{B^)^r]ri 


63.5 


58.4 


63.7 


58.2 


51.4 


46.5 


B^{B^)^r]'r]' 


69.9 


78.4 


74.0 


77.1 


76.9 


62.9 


if {B^)^rir]' 


71.0 


66.7 


75.5 


73.2 


67.3 


69.2 


B\B^)^K^7r^ 


-6.53 


-16.1 


-5.91 


-14.8 


-4.20 


-10.7 




37.8 


38.0 


36.7 


36.9 


29.3 


29.4 




-7.24 


-15.0 


-6.69 


-14.1 


-4.83 


-10.3 




-0.84 


-1.05 


-0.99 


-1.24 


-0.88 


-1.10 


B^^K^T] 


8.35 


14.0 


12.8 


21.1 


16.3 


26.3 


B'^^K^T]' 


-3.19 


-7.50 


-3.41 


-8.12 


-2.75 


-6.62 


B^{B^)^K'^r] 


39.1 


39.2 


38.2 


38.4 


30.7 


30.8 


-^(B^)^K^r]' 


34.0 


33.9 


32.0 


32.0 


24.8 


24.8 


B^^K^K'^S 


21.0 


26.4 


16.6 


20.2 


10.9 


13.1 




18.0 


15.4 


13.8 


11.5 


8.97 


7.38 


'^{B^)^K+K- 




34.7 




1.65 




-32.4 
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TABLE IX. The dependence of CP-averaged branching ratio and CP asymmetries for 
B^^K^TT^ on Q, (p in endpoint divergence fl dx/x. BRf and BR^+'^ have the same mean- 
ings as those in Table IV. A^p and can be defined in a similar way. The results are 



calculated with default values of other parameters at /i = nib using the BENS approach. 



Q 







3 










6 




0(deg) 


0~360 





90 


180 


270 





90 


180 


270 


BR-^(10-6) 


6.37 


6.31 


6.36 


6.43 


6.38 


6.25 


6.35 


6.49 


6.39 


BR^+"(10-'5) 


7.51 


11.8 


6.75 


6.43 


6.18 


21.3 


4.63 


7.89 


3.48 




8.28 


8.45 


7.15 


8.11 


9.40 


8.63 


6.02 


7.94 


10.5 




7.24 


5.15 


14.2 


8.11 


1.47 


3.32 


29.5 


6.81 


-12.4 



30 



FIGURES 




(a) (b) (c) (d) 




(e) (f) (g) (h) 



FIG. 1. Order of ag corrections to hard-scattering kernels. The upward quark hnes represent 
the ejected Hght meson from b quark weak decays. These diagrams are commonly called vertex 
corrections, penguin corrections and hard spectator diagrams for Fig.(a)-(d), (e)-(f) and (g)-(h), 
respectively. 




FIG. 2. Order of Og corrections to annihilation diagrams for B^PP 
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FIG. 3. CP-averaged branching ratios as functions of 7. The soUd lines are drawn with 
the default values at the scale of /i = rrib, the horizontal slashed-line bands correspond to 
experimental data within one standard error, and the dot-shades denote the variation of the 
theory input parameters, including the CKM elements, the form factors, the uncertainty from 
weak annihilation and hard spectator scattering. 
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FIG. 4. The dependence of CP-averaged branching ratios for ir^K^, B^^K^ir^ 

and vr+vr^ on the variations of input parameters under the QCDF approach at the 

renormaHzation scale /i = m-f,. The hues and bands have the same meaning as in Fig. 3, and 
the dot-shades correspond to the variation of the CKM matrix elements, in particular \Vub\ (the 
first row), form factors (the second row), endpoint divergences in the hard spectator scattering 

(the third row) and weak annihilations X^ithe forth row). 
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FIG. 5. Ratios of CP-averaged branching fractions versus 7. The hnes, bands and dot-shades 
have the same meaning as in Fig. 3. 
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